Magnetic Phase Transition in FeRh 
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Abstract 

Density functional calculations are performed to investigate the phase transition in FeRh alloy. The effective exchange 
coupling, the critical temperature of magnetic phase transition and the adiabatic spin wave spectrum have been obtained. 
Different contributions to the free energy of different phases are estimated, ft has been found that the antiferro-ferromagnetic 
transition in FeRh occurs mostly due to the spin wave excitations. 



The antiferromagnetic (AFM) - ferromagnetic (FM) 
phase transition in FeRh with the ordered CsCl struc- 
ture has been intensively studied both experimentally 
and theoretically. This transition occurs at Ttr ~ 340 
K without any accompanying structural changes 0, al- 
though there is an abrupt 1% volume expansion. At low 
temperatures the magnetic configuration of FeRh is the 
type-II AFM (successive layers of (111) Fe planes AFM 
coupled) with moments 3.3/LtB on Fe atoms (Rh atoms 
are nonmagnetic) . Above Ttr in FM phase magnetic mo- 
ments are 3.2/iB on Fe atoms, l.O/is on Rh atoms and 
the Curie temperature (Tc) is « 670 K [||. In addition, 
it was found that Ttr is increased with pressure 0, 0|. 
The AFM-FM transition can also be induced by apply- 
ing external magnetic field, whose critical value at zero 
temperature is about 300 kOe and smaller at higher tem- 
peratures 0, making this material a natural magnetic 
multilayer with a large magnetoresistance effect. 

Early theories of this transition based on the exchange- 
inversion model 6] , which assumes a change of sign for 
the exchange parameter at some volume, can not ac- 
count for the experimental observation of the large en- 
tropy changes at Ttr |Ei IE Hi 1^ - After discovering that 
the electronic specific heat in the FM phase is ne arly 
four times larger than in the AFM phase, Tu et al. [13 
proposed that the change of the band electron entropy 
plays a major role in this transition. However, they used 
iron-rich alloys where the value of the specific heat is very 
sensitive to the concentration jllj] . Consecutive dielectric 
function measurements demonstrated that the band 
structure of FeRh is not drastically modified during the 
AFM-FM phase transition. 

First-principles band-structure calculations were also 
carried out to study this transition. Earlier calculations 
[Tsf did not compare the relative stability of the AFM 
and FM states. Moruzzi and Marcus [3| confirmed that 
the type-II AFM structure is the ground state, while the 
FM structure represents another stable solution with the 
total energy nearly /S.E=2 mRy/atom higher at a larger 
volume. Similar results were obtained in Refs.^Tsl Il6j|. 
However, the energy difference IS.E obtained in these cal- 
culations appears to be much larger than the experimen- 
tal data and, overall these studies did not provide any 
convincing explanation of the nature of this transition. 
The experimental AS, deduced either from the latent 
heat TtrAS* or from the critical magnetic field at zero 



temperature Q, is about 0.2 mRy/atom, which is an or- 
der of magnitude smaller than the calculated value. On 
the other hand, while in Ref.^j the calculated equilib- 
rium lattice constant of the FM state is only 0.5% larger 
than in AFM one, the energy of latter remains lower than 
the energy of the FM state until the lattice constant is 
increased by 3%. The authors 14] proposed that the zero- 
point lattice vibrations can correct the total energy re- 
sult. However, from the Debye temperature Qd calcu- 
lated in that paper, one can find that the correction due 
to zero-temperature vibration energy = dkBOo/^ is 
nearly two orders of magnitude smaller than the calcu- 
lated AE. Gruner et al. investigated the thermo- 
dynamic behavior of the system within Ising model and 
found that at Ttr the free energy, gained due to thermoex- 
citation in the FM state, is nearly 0.02 mRy/atom larger 
than that in the AFM state. So, they proposed that it 
is this thermocxcitation that drives the transition. How- 
ever, the magnitude of their free energy change appears 
too small to compensate the internal energy loss for the 
transition to occur. So far, in spite of many years of re- 
search, there has been no convincing explanation of the 
nature of the phase transition in FeRh. 

In this paper we study this transition using first prin- 
ciples calculations and the non-collinear version of the 
linear muffin-tin orbital method in the atomic-sphere ap- 
proximation (LMTO-ASA). In the local spin density ap- 
proximation (LSDA) we use the von Barth-Hedin po- 
tential, and for the nonlocal corrections the Langreth- 
Mehl-Hu functional ,2&|, the relativistic effects and the 
combined corrections |l3 are included. For the radii of 
atoms, the ratio R^h/Rpe = 1-03 was used. The self- 
consistent calculations are performed for the different lat- 
tice parameters for FM and AFM states using a spin spi- 
ral approach. The gradient corrections are expected to 
be important in Fe-rich BCC based systems due to well 
known fact that LSDA predicts, for instance, the wrong 
FCC ground state pure Fein. Below we will show that 
it is exactly a case in FeRh. Using electronic density of 
states (DOS) and exchange parameters, we calculate the 
free energy change of the AFM and FM states. Vari- 
ous thermal quantities related to the AFM-FM transi- 
tion, including the transition temperature, its pressure 
dependence, the entropy and the specific heat changes 
are calculated and compared with the corresponding ex- 
periments. According to our calculations, the AFM-FM 
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transition in FeRh appears primarily due to the magnon 
(spin wave) excitations. 

Fig. 1 (a) shows the calculated energy of the AFM 
and FM states. In LSDA the equilibrium Wigner-Seitz 
radius i?ws for the AFM (FM) phase is 2.767 (2.780) 
a.u., which is smaller than the previous non-relativistic 
result 2.782 (2.798) a.u.Q| (our non-relativistic result is 
2.789 (2.803) a.u.). The energy difference between the 
AFM and FM states at their respective equilibrium Rws 
is 1.89 mRy/atom. With nonlocal corrections this en- 
ergy difference is reduced to 0.206 mRy/atom, which is in 
agreement with the experimental value 0.196 mRy/atom 
13 • Other results are listed in Table I. 

Only coUinear AFM and FM configurations were ob- 
served in the experiments. In our calculation, however, 
we found that without nonlocal corrections the FM con- 
figuration is not a locally stable state with respect to 
the magnetic moment deviations, as is shown in Fig. 1 
(b) (two kinds of noncoUinear configurations are consid- 
ered) . In both those states the Fe atoms are divided into 
two sublattices, identical to that in the type-I (successive 
(001) Fe layers belong to the different sublattices) and 
type-II AFM states with the angle 9 between Fe mag- 
netic moments from the different sublattices, while the 
moments of Rh atoms are parallel to the sum of the Fe 
moments. Such instability of the FM state is removed 
if the nonlocal corrections are taken into account. Re- 
sults from Fig. 1 (a) and (b) suggest that the nonlocal 
corrections are important in FeRh and necessarily should 
be taken into account in the investigation of the phase 
transition. 

Using formalism from Ref.|^ we calculated the pa- 
rameters of the exchange coupling. The obtained 
parameters are different for the FM and AFM con- 
figurations (see Table II). To estimate the mag- 
netic contributions to the free energy we calcu- 
lated the magnon spectrum in both phases. In 
the adiabatic approximation for the AFM state 

^AFM ^ (2g/iB/mFo)^(Jo - Jq)(Jo " Jci+q), with Jq = 

12 j Jije^'^'^'^ being the Fourier transformation of Jy 
in AFM state and Q = (tttttt). For the FM state 
w™'± = 5MB (^q + Sq ± y^(Aq-Bq)2+4X2), whcrc 

= (JF^F^ - J^<=F<= + 4<^^^)/mFe, = (J«h Rh „ 

^ jFcRh)/„jj^j^^ ^j^d = Jf^VV"^FcmRh, with 

jFcFc^ jFcRh^ g^^^ j^iiT^h being the Fourier transforma- 
tions of the exchange interactions inside (or between) the 
corresponding sublattice(s). The magnon DOS obtained 
from this spectrum is shown in Fig. 2 together with the 
electronic DOS. 

To study the relative stability of the AFM and FM 
configurations at finite temperatures, let us compare 
their free energies. Here we consider the contributions 
from electrons and magnons only. The lattice contri- 
bution is neglected because the magnitudes of the bulk 
moduli and Debye temperatures in both phases are very 
similar (Table I). The free energy due to band electrons 



and magnons are given by 
FeiiT) - ^{epn - keT J deN{e) log(l + e^^^-^^/'^'^^)} , 

WT) = -^/(03Ml-e--/^-^), (1) 

where n and N{e) are the number of electrons and the 
electronic DOS, correspondingly. Fig. 3 (a) shows the 
free energy difference AF = — F^^^ as a func- 

tion of temperature. The transition temperature Ttn 
determined from AF — 0, is 371 K, which is close to 
the experimental result Ttr ~ 340 K. Both contributions 
are shown, with the main one (more than 80%) coming 
from magnons. So, the origin of the AFM-FM transition 
should be attributed primarily to the magnon excitations 
rather than to pure electronic spectrum modifications. 
The obtained AF{T) also enables us to get the pres- 
sure dependence of the transition temperature, whose ex- 
perimental value is about dTtr/dP « 5.1 ~ 5.8 K/kbar 
[3, 01 . With applied pressure the AFM state gains more 
free energy AG = PAV {AV is the volume difference) 
than the FM state and the derivative dAG/dP is close 
to —6.2 X 10^'^mRy/kbar per atom. From Fig. 3 (a) 
dAF/dT = -1.17xl0-3mRy/K at Ttr and from the equi- 
librium condition dAF = dAG we obtain dTti/dP w 5.3 
K/kbar, which agrees well with the experiments. 

From the obtained electronic DOS and magnon spec- 
trum, one can also evaluate various thermal quantities. 
In Fig. 3 (b) we show the calculated differences of the 
entropy and the specific heat between AFM and FM 
states as a function of temperature. These two quanti- 
ties are independent of the zero temperature energy, with 
their measured values at Ttr being AS^p « 13 — 19.6 
Jkg-iK-i|l|il3 and AC'^P w 13.6 Jkg-^K-i 0. Just 
as the free energy, the calculated contributions to both 
AS* and AC near Ttr mainly determined by the magnon 
excitations. At the calculated Ttr — 371 K, we obtain 
AS" = 19.3 Jkg-iK^i and AG = 15.1 3kg-^K-\ while 
at Ttr ~ 340 K, the corresponding values are 17.9 and 
15.6 Jkg-iR-i. 

Since the Rh atoms have a nonzero value (l/is) of mag- 
netic moments in the FM state, this state has more mag- 
netic degrees of freedom. It was proposed that this addi- 
tional number of degrees of freedom increases the entropy 
and thus stabilizes the FM state 0|. The evaluation of 
this entropy gain gives AS* w Nks log 2 « 36 Jkg~^K~^, 
which is too large compared to the experimental results. 
According to our calculated DOS of magnon, this picture 
of extra entropy is not quite accurate. /^From Fig. 2 (b) 
one can see that in the FM state the magnon excitations 
cjq and are separated by the large energy gap. We 
call the cjq i^q) mode to be Fe(Rh)-like, because a simi- 
lar mode (the dotted line in Fig. 2(b)) can be obtained if 
we fix the orientations of the Rh (Fe) moments by letting 
Xq = in the expression for uj™. Near Ttr only the Fe- 
like mode contributes to the thermal properties so that 
the number of the effective magnon states in the FM and 
AFM phases is the same, i.e., there is essentially the same 
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number of spin degrees of freedom in the FM and AFM 
states. That does not mean, however, that the nonzero 
Rh moments in the FM phase do not contribute to the 
thermal properties. On Fig. 2(b) it is shown that with- 
out the movement of the Rh moments there is the energy 
gap between the ground state and the lowest-energy Fe- 
like magnons uj„—q . A calculation shows that at Tti the 
magnon free energy (with the orientations of the Rh mo- 
ments being fixed) is only about one third of that when 
they are not fixed. In other words, the Rh moments 
soften considerably the stiffness of the Fe-like magnons, 
significantly influencing the thermal properties. A com- 
parison of the magnon DOS in the AFM and FM states 
indicates that near Tn it is much easier to excite the Fe- 
like magnons in the FM state. This is the main reason 
for the difference in thermal properties between the two 
phases, and it is also the driving force of the AFM-FM 
transition in FeRh. 

Finally let us evaluate the Curie temperature by using 
the obtained pair exchange interactions Jij . In the mean 
field (MF) approximation 
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From Jo*«*« = -3.56 mRy, Jg^^^*^ = 9.92 mRy and 
jRhRh ^ 0.85 mRy we obtain 927K, which is nearly 40% 
higher than the experimental Tc ^ 670K. Our Monte- 
Carlo calculations (with all calculated long ranged Jij 
included) produced correspondingly 660 — 690 K, so that 
the ratio Tc/T^^ w 0.71 - 0.74 is nearly the same as 
that of the simple cubic lattice Heisenberg model where 
Tc/T^^ = 0.722 J3J. The agreement of the calculated 
and the experimental Tc indicates that the Heisenberg 
model may still work well in the temperature region near 
Tc in FeRh. 



FcFc I rRhRh 
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Fig.l. The total energy obtained in the local (L) and 

nonlocal (NL) approximations, represented by open and 

solid symbols. The total energy for: (a) AFM (circles) 

and FM (squares) states as a function of i?ws; (b) for the 

type-I (squares) and type-II (circles) noncollinear states 

as a function of the spin spiral angle 6. Rws is fixed at 

the equilibrium value of the FM state. 

Fig.2. Calculated DOS of FM (solid line) and AFM 

(dashed line) states: (a) the electronic DOS and (b) the 

magnon DOS. In (b) the dotted lines correspond to the 

FM magnon DOS when X„ = in the expression for 

w™. 
q 

Fig. 3. The calculated differences of (a) the free en- 
ergy AF, (b) the entropy AS* and the specific heat AC 
between FM and AFM phases. The dashed and dot- 
ted lines (short-dashed and short-dotted lines for AC) 
correspond to the contributions from the electrons and 
magnons, correspondingly. The solid line is their sum. 



TABLE II: Pair exchange parameters (in mRy) in AFM and 
FM phases of FeRh. Corresponding coordinates are shown in 
units of the lattice constant. 



Type of pair Aa; Ay Az 


tFM 

<iij 




Fe-Rh 


0.5 


0.5 


0.5 


1.062 







1.5 


0.5 


0.5 


0.058 







1 








-0.098 0.442 


Fe-Fe 


1 


1 





0.104 


0.008 




1 


1 


1 


-0.479 


0.603 




2 








0.120 


0.099 




2 


1 





0.045 


0.005 


Rh-Rh 


1 








0.086 







1 


1 





0.018 






TABLE I: Calculated physical properties of the AFM and 
FM configurations of FeRh obtained in the local (first row) 

and nonlocal approximations, mpo and rriRh arc magnetic 
moments of Fc and Rh atoms, B is the bulk modulus, 0d is 
the Debye temperature and N{tF) is DOS per formula unit 
at the Fermi level. 

Tiws AE mpo niRh B 9o Niep) 

(a.u.) (mRy/atom) (mb) (M-b) (kbar) (K) (Statcs/Ry) 

2.767(2.780) 1.89 3.12(3.22) 0(1.04) 2454(2364) 385(379) 18.0(29.5) 

2.796(2.807) 0.206 3.28(3.31) 0(1.02) 2194(2181) 366(365) 15.6(28.0) 
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